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Abstract 

We introduce an equation named matrix Dirac equation which can 
be considered as a generalization of Dirac equation for an electron. The 
liaison between matrix Dirac equation and standard Dirac equation is 
discussed. We write a lagrangian from which matrix Dirac equation can be 
derived. This lagrangian is invariant under global unitary transformations 
of variables. The requirement of a local (gauge) invariance of lagrangian 
leads us to lagrangian with gauge fields. 

Introduction. 

After the famous Dirac equation for an electron was found in 1928 |§, works 
of H.Weil Q, V.A.Fock and others appeared, in which electromagnetic field 
was described as gauge field of the Dirac equation appearing from the demand 
of local gauge invariance with respect to phase transformation (an abelian U(l) 
gauge group ) of Dirac's lagrangian. Further development of this approach for 
the non-Abelian gauge fields was made in the work of Yang and Mills (1954) Q 
considering the group of isotopic transformations (gauge group SU(2)). Their 
work was soon generalized for the wide class of Lie groups. Non-Abelian gauge 
fields began to be named Yang-Mills fields, whereas the equations describing 
them were named Yang-Mills equations. In modern physics such fields are used 
in models of electroweak and strong interactions. 

In the present article an equation named matrix Dirac equation is intro- 
duced. This equation can be considered as a generalization of Dirac equation 
for an electron. Certain features of matrix Dirac equation are investigated (cur- 
rents, canonical forms). Further on the basis of local gauge invariance regarding 
unitary group a system of equations is introduced consisting of matrix Dirac 
equation and equations of Yang-Mills or Maxwell. This system of equations de- 
scribes Dirac's field interacting with the gauge field of Yang-Mills or Maxwell. 
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Some ideas proposed in [|[ ^, |l^ are used in the article. Certain elements of 
the construction proposed could be found in the works of Hestenes ^ , Kahler 
[|, Pestov 111, Pezzagha and Differ [[o|. 

1 A standard Dirac equation. 

Let us consider a Dirac equation and its generalizations. A vector x = {x'^ , x^ , x^, x^) € 
72.^ defines a point in space-time, x'^ - time coordinate, x^,x'^,x^ - space coor- 
dinates, = d/dx'^, /i = 0, 1, 2, 3 are the partial derivatives and □ = d^d^ = 
do^ — di^ — d'l' — dz' . A Klein-Gordon equation for a function (\) — 4>{x) 

{D + m^)(t) = 0, (1) 

where m - nonnegative real number, describes particles with spin and mass 
m. For a description of spin 1/2 particles (electrons) P. A. M. Dirac suggested 
a system of equations of first order that he received by factorizing the Klein- 
Gordon operator 

{i-f^'df,+m){i-f''d^ -m) ^ -{a + 771^), (2) 
where 7^, fj, — 0,1,2,3 are algebraic objects that satisfy the relations 

7''7'^+7V-2g^^ ^,^ = 0,1,2,3 (3) 

with Minkowski tensor g — {g^'^) — diag(l, —1, —1, —1). 7'' can be represented 
by the matrices of order no less than four. We shall use the following represen- 
tation for 7^: 

which is called the Dirac representation [a^ - Pauli matrices). Let 1 denote the 
identity 4x4-matrix and 7^^^ = 7^7"^ for < ^ < < 3, 7^'''^ = ^t^^^j^ for 
0</i<i/<A<3, 7^== 7"i23 = -/°-/'^-/'^-/^. The 16 matrices 

1, 7^ r", f (5) 

are linear independent and form a basis of A^(4, C) - algebra of four dimensional 
complex matrices. 

The Dirac equation has a form 

(17^9^ - ml)V' = 0, (6) 
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and if we consider the aggregate in brackets as a 4x4-inatrix, then tp — ip^x) 
must be a matrix with four lines and an arbitrary number of columns. The 
identity (^) asserts that all the components of the matrix ip (if they are smooth 
enough) satisfy the Klein-Gordon equation. 

We shall call the equation (j^) with one column matrix the standard Dirac 
equation. In that case ^ is called bispinor or Dirac spinor. One can also consider 
the equations (^) with / > 1 columns in the matrix tp. For the different purposes 
physicists use the Dirac equation with the different numbers of columns in 
ip. The following list does not lay claim for completeness or indisputability: 

/ = 1: Quantum electrodynamics (a gauge field theory with U(l) symmetry 
group). 

/ = 2: Theory of electroweak interactions (SU(2) gauge field theory). 

I — 3: Theory of strong interactions - quantum chromodynamics ( SU(3) gauge 
field theory). 

I = 4: The Dirac equation with ^ = 4 is called a Rarita-Schwinger equation. 
It is used for a description of spin 3/2 particles. 

^ = 5: Georgi and Glashow have suggested SU(5) gauge field theory as 
Grand Unified Theory (GUT). 

I > 6: SU(Z) gauge field theories are developed by physicists as GUT. 

Gelfand, Minlos and Shapiro considered all relativistic invariant systems of 
equations of the form 

L''^^^P + ini) = 0, (7) 

where are square matrices, i/) - vector, k - real constant. 
There is an evident generalization of the identity (^) 

+ m(zl - 2/7'^))(z7"a, - m(zl + y-f')) = -(□ + m^), (8) 

where z, y are complex constants and + = 1. It leads to the equation 

(z7^9^-m(zl + y75))7/> = 0, (9) 

We can consider the factorizations (^,(||) of the Klein-Gordon operator as one 
of possible methods of a reduction of the Klein-Gordon equation of second order 
to a system of equations of first order. There is another method of such a 
reduction that leads to the following system of equations of first order: 

- m{iPN + j^^K) = 0, (10) 

which depends on two matrices N,K ^ Ai{l,C). We shall call the equation 
(lOf) the matrix Dirac equation, emphasizing that unknown (wave) function 
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tp = ip{x) in ( [L0[ ) is a 4xZ-niatrix. If matrix ip has only one column, then we 
get an equation (^), or standard Dirac equation (when N = 1,K = 0). The 
equation (p^), generally speaking, can't be reduced to the equation of the form 

(0)- 

Theorem 1 If a Ax l-matrix ip — tpix) with twice continuously differentiable 
elements is a solution of where the matrices N,K G ^A{l, C) do not depend 
on X and satisfy the relations 

[N,K] = NK - KN ^0, N^+K^ = li, (11) 

where 1; - is identity I x l-matrix, then the matrix ijj is also a solution of the 
Klein- Gordon equation. 

Proof. Let us consider an action of the operator ij'^d^ from left to the equation 
® 

{i-f'^df.fil; - m{ijt'df,'ip)N + m-f^{ij>'df,'iJj)K = 0, 



and use the relations (|11|) and 

(i7^a^)2 ^ i-f'^df.'ip = m{iPN + ^i^iPK). 

As a result we get the Klein-Gordon equation (□ + m?)'4> — 0. Theorem is 
proved. 

The formula ( p^ gives us a set of equations that depend on two matrices 
iV, K with the relations (pi]). How to describe all matrix pairs iV, K that satisfy 
(pi])? If we write the matrix N using a Jordan normal form, then we can find 
all the corresponding matrices K by doing an elementary calculation |l^. For 
example, when I = A we get 15 classes of pairs N,K that depend on several 
parameters. 

Example 1. The matrices 

iV = Fdiag(2i,...,zO^"\ K^VA:mg{y^,...,yi)V-\ (12) 
where z^^ + yk^ = 1, V- nondegenerate matrix from A4{l,C), satisfy (pi]). 
Example 2. Let 1 = 4. The matrices 

/ y a b c ^ 




^ = ^ n n ; V-\ K = V 



y a b 
y a 
VO 



v-\ 



where + = 1, y ^ 0, a = -z/y, b = -l/{2y^), c = -z/{2y^) and V- 
nondegenerate matrix from A4{A,C), satisfy (plf). 
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Let us consider a solution ip of (|l^) with the matrices TV, K from the example 

1 and connect it with the solutions of the standard Dirac equation. The columns 
of the matrix V denote hy Vk {k — 1, . . . ,1). Then Vk are the eigenvectors of N 
corresponding to the eigenvalues Zk- Simultaneously, Vk are the eigenvectors of 
K corresponding to the eigenvalues tk- So, if we multiply from right by Vk 
and denote Tpk = V'^fc then we come to I equations 

iYd^ipk - m{zkl + ykl^)i^k = 0, fc = 1, . . . , / 

that have the form or (^). 

2 A lagrangian of the matrix Dirac equation. 

The standard Dirac equation (|^) can be derived from Dirac's lagrangian 

(i/2)(Vi7^a^^ - 5,,VW') - m^V- (13) 

with the aid of variational principle . N,K are hermitian matrices from 
M.{l^C), then the matrix Dirac equation ( p^ can be derived with the aid of 
variational principle from the lagrangian 

(l/2)tr (i(V'7^9^?/' - 9^V'7''V') - m{^^N + N^P^) + i^-^^i^K + K^!-/^^)) , 

(14) 

in which tr is a trace of matrix and ip = ip'^j'^. 

Let us introduce a commutator algebra of matrices N, K: 

com{N, K)^{V e Mil, C) : [V, N] = [V, K] = 0} 

and a group 

g{l, N, K) = coni(Ar, K) n U(0, 

where U(Z) is a group of unitary matrices from M.{l,C). The group Q{l,N, K) is 
a compact Lie group. A set of all antihcrmitian matrices (V^^ = ~V) commuting 
with N and K 

C{1, N, K) = com(iV, K) n u(0. 

can be considered as a real Lie algebra of the Lie group Q{l,N, K). 

The lagrangian is invariant under global (not dependent on x) transfor- 
mations 

^^ipV, V-^ip, with V e 0(1, N, K), (15) 

This fact is evident from the identity ti-{V^^BV) — trB. In order to get a 
lagrangian which is invariant under a local (gauge) transformations ( [l5| ) , where 
V = V{x) is a function of x with values in Q{l,N,K), we must replace in 
( p4[ ) partial derivatives dfj,ip, d^i^ by the respective covariant derivatives D^ip = 
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d^ifj — ipttfi,, D^if) = dfj_ip + that depend on functions = a^(a;) with its 
values in the Lie algebra C{1, N, K). A transformation rule for has a form 

a^,^V-^a^,V + V-^d^V, Veg{l,N,K). (16) 

There is a complete gauge invariant lagrangian, which describes a field "0 
interacting with the gauge field a^, 

L = ^tT{i{j;-/^'{^|,iJ - -iAa^J - (9^-0 + 0^0)7^1/') 7^ 
-m{i:iPN + NtPi: + ^Pj^i^K + K^Pj^^P)) + itr(/^^/^''), ^ ' 

where /^j/ = d^a^, — d^a^i — [ajy, a^,]. If N — N'' , K — K\ then the lagrangian 
L leads to the system of equations 

h''{df,ip - V'fl^) - m{^PN + j^^PK) = 0, 

d^a„ - d^a^ - [a^,a^] - f^^ = 0, (ig) 

c 

7r: u(Z) jC{1,N,K) is a projector operator to the Lie algebra C — C{l,N,K) 
which can be considered as a linear subspacc of the vector space u(^) of an- 
tihcrmitian matrices. The second and the third equations in (18) are called 
Yang-Mills equations. 

Definition. If y = V{x) is a function of x with the value in a Lie group 
g C Q{1, N, K), then the transformation 

^ a'^ = V-'a^V + V-^d^V, (19) 

is called a gauge transformation of fields ip, a^, f^u with the group Q (note, that 
the third relation from (|l9|) is a consequence of the second). 



The lagrangian ( |17D and the system of equations ( |18| ) are invariant under a 
gauge transformation (|l^). That means L{tp^a^) = L{tp',a'^), and if tp^a^jf^i, 
satisfy then ip'^a'^J^^ from (|l|) also satisfy (|l|). 

Theorem 2 Let N,K e M{l,C) be such, that 

= a2l(+/32P2, P2Pj-l;, Pj=P2, a2,/?2e7^, 



and = ip{x) is a solution of the matrix Dirac equation Let us denote 

J" =^ {i^Y^), where C = C{l,N,K) = com{N , K) (1 u{l) . Then 

d,r{iP) = 0. (21) 



6 



In other words, if we pose additional conditions ( pO| ) on N,K, then a right 
hand part of Yang-Mills equations is a current of matrix Dirac equation. 

Proof. A solution ip of the matrix Dirac equation ( p^ also satisfies an identity 

5^(iV^7^V) - mii^^N - N^il) + 7/57^ - K^-f^i:) = 0. (22) 

which is consequence of (p^). To show this we must multiply ([To|) from left 
on 'ijj and subtract hermitian conjugated equation — i9^V^(7^)^ ~ m{N^ifj^ — 
K''' tjj'' j^) = 0, multiplied from right on j^tp. The result can be written in the 
form (|2|). 

Let us prove that if K satisfy (|2^) , then 

n {■4>%l)N - + ij-f^ijK - K^-f^ij)) = 0. (23) 

It is easy to check, that the matrices 

anticommute with the matrices Pi , P2 respectively: 

BiPi = -PiBi, B2P2 = -P2B2. 

We can rewrite this fact as 

Bi e M{l,C)\com{Pi), B2eM{l,C)\com{P2), (24) 

where com(P) C A4{l,C) is a subspace of all matrices commuting with the 
matrix P e M{l,C). It follows from (|2|), that 

B1+B2 e7W(/,C)\(com(Fi)ncom(P2)). (25) 

Taking into account the definition of the Lie algebra C, we get from ( p5| ) the 

C c 
relation tt (Bi + B2) — 0. Acting by the projector operator tt on the right hand 

and left hand parts of the identity (^) , and using the formula ( p3| ) , we get 

^(9^(z^7^V)) ^S^J^^O. 

Theorem is proved. 



3 A general form of the matrices A^, K. 

The matrix Dirac equation (^) depends on two matrices N, K. Let us compile 
all conditions on the matrices K. 
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If N, K satisfy (^T|), then the solution -0 of equation (|l^) also satisfies Klein- 
Gordon equation. 

If iV, K are hermitian, then the matrix Dirac equation can be derived from 
the lagrangian ( 14 ) with the aid of variational principle. 

If iV, K satisfy ( ^ ) , then the right hand part of Yang-Mills equations 
is a current of matrix Dirac equation. 

Evidently, if N, K satisfy (|2^), then K are hermitian matrices. 



Theorem 3 The matrices N,K e M{l,C) satisfy ([Hp and ^2(\ ) if and only if 
N — C/tdiag(cos^, . . . ,cos^,cos77, . . . ,cos77)J7, 



p pieces q pieces 

K = J7'''diag(sin^, . . . ,sin^,sinr/, . . . ,sin?7)[/, 

v ' " V ' 

p pieces q pieces 



or 



N ^ C/tdiag(±l,...,±l)C/cosC, 
K = C/tdiag(±l,...,±l)C/sin^, 

where p + q — I, U -unitary matrix and < ^, < 27r. 



(26) 



(27) 



Proof. The matrix N in ( pO| ) depends on a hermitian and simultaneously 
unitary matrix Pi. So, Pi can be written with the aid of a unitary matrix U: 

Pi = f7tdiag(±l,...,±l);7. 

Hence 

iV = [/tdiag(Ai,...,AOC/, 

where = A+ = ai -I- /3i, or Xk = = ai — Pi. The same is true for the 
matrix K: 

K = V^diag{ei,...,ei)V, 

where V^V = 1, = e+ = q;2 + 02, or ek = e^ ~ a2 — ^2- Matrices N,K 
commute with each other and can be simultaneously reduced to a diagonal 
form with the same similarity transformation jl^ . That means, that we can 
take U = V. The condition N'^ + K"^ = 1 gives Xk'^ + Sk'^ = I, fc = 1, . . . , I. 
There are two possibilities. The first: (A+)^ ^ (-^ i^^T )^ and so 

Afe = A+ = cos^, efc = e+ = sin^, k ^ 1,. . . ,p, 

Xk = X^ = cosi], ek = = sinrj, k = p + 1, . . . ,1. 
The second: (A+)2 = (A-)^, (e+)2 = (e-)2 and so 

Afc = ±cos^, efc = ±sin^, fc=l,...,L 

Theorem is proved. 
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4 A polar gauge. 



There is a theorem about the polar decomposition of matrix. 

Theorem 4 ([^^) An arbitrary matrix M G jVl(l^C) can be written in a form 
M — PU , where P £ A4{l,C) is a hermitian nonnegative defined matrix with 
the same rank as M , and a matrix U G M.{l^C) is unitary. 

Let us consider a system of equations ( |l8| ) where 1 — 4:, N — Icos^, K = 
Isin^ with the gauge group G = U(4). And let ij^^a^jf^^i, be a continuously 
differentiable functions of x which satisfy and such that ij) G A4(4,C) and 
Opi/^ij/ G £ = u(4). Using the theorem 4, we get that in every point x G TZ'^ 
the matrix -0 = ip{x) can be written in a form ip = PU, where P E A4{1, C) is a 
hermitian nonnegative defined matrix and [/ G A4(/, C) is a unitary matrix. 

Let us suppose, that the matrix ^ = 'ip{x) can be written in a form il){x) = 
P{x)U{x) in some region il C TZ^ and P{x),U{x) have continuously differen- 
tiable elements for all x G fi. In that case the solution ip = ip{x) of (|8|) in the 
region fl defines a unitary matrix U = U (x) with continuously differentiable ele- 
ments. We can take V — V{x) — U~^. The system of equations ( p^ ) is invariant 
under a gauge group U(4), and after a gauge transformation with V G U(4), we 
come to the solution of ( p^ ) 

^' = ipV. a'^ = V-^a,y + V-^d^,V, F' = V^FV, 

where rp' = ip'ix) is a hermitian nonnegative defined matrix and we name it 
solution of (|l8|) in polar gauge. 

In quantum mechanics particles are described by the wave functions which 
belong to some complex finite or infinite dimensional vector (hilbert) space. 
Observables are hermitian operators on that space. It will be naturally to 
suppose, that the solution of (|l8|) is a wave function of fermion. 

If we take ip in the polar gauge (in the case of gauge group U(4)), then ip is 
a hermitian 4x4-matrix equivalent to the hermitian 16 x 16- matrix ^' (X) 1. That 
means, that we can consider ip not only as a wave function, but simultaneously 
as some observable value. 

Let us make two final remarks. 

It can be shown, that the matrix Dirac equation (^|) is invariant under the 
Lorenz transformations. 

In case / = 4, the wave function ip, which satisfy ([ic|), is a 4 x 4-matrix 
and so, it can be represented as a linear combination of basis vectors (^) of 
Dirac's algebra (Clifford's algebra). Substituting it into ([l0|), we get an equation 
written in Clifford's algebra terms and not dependent on matrix representations. 
Such an equation can be naturally generalized to the curved space-time with an 
arbitrary pseudorimanian metric. 

This questions will be considered in the next publications. 
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